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ABSTRACT 

It is generally believed that O stars, confined near the galactic midplane, are somehow able to 
photoionize a significant fi-action of what is termed the "diffuse ionized gas" (DIG) of spiral galaxies, 
which can extend up to l-2kpc above the galactic midplane. The heating of the DIG remains 
poorly understood, however, as simple photoionization models do not reproduce the observed line 
ratio correlations well or the DIG temperature. We present turbulent mixing layer models in which 
warm photoionized condensations are immersed in a hot supersonic wind. Turbulent dissipation and 
mixing generate an intermediate region where the gas is accelerated, heated and mixed. The emission 
spectrum of such layers are compared with observations of Rand (ApJ 462, 712) of the DIG in the 
edge-on spiral NGC891. We generate two sequence of models that fit the line ratio correlations 
between [Sii]/Ha, [OiJ/Ha, [Nii]/[Sii] and [Oiii]/II/3 reasonably well. In one sequence of models the 
hot wind velocity increases while in the other the ionization parameter and layer opacity increases. 
Despite the success of the mixing layer models, the overall efficiency in reprocessing the stellar UV is 
much too low, much less than 1%, which compels us to reject the TML model in its present form. 
Subject headings: Line: formation - Hydrodynamics ~ Turbulence - Galaxies: individual: NGC 891- 
ISM: lines and bands 



1. INTRODUCTION 

The vast majority of the free electrons in the ISM of 
the Milky Way reside in a thick 900 pc scale height) 
diffuse layer, which fills about 20% of the ISM volume, 
with a local midplane density of about 0.1 cm"'^. Such 
a phase is now known to be a general feature of ex- 
ternal star- forming galaxies, both spirals (e.g., Walter- 
bos 1998; Otte etal. 2001, 2002) and irregulars (e.g.. 
Hunter & Gallagher 1990; Martin 1997). It is referred 
to as the warm ionized medium (WIM) or the diffuse 
ionized gas (DIG). For edge-on spirals, only in the more 
actively star-forming galaxies does the gas manifest it- 
self as a smooth, widespread layer of emission detectable 
above the Hii region layer (Rand 1996). It is generally 
believed that O stars, confined primarily to widely sep- 
arated stellar associations near the Galactic midplane, 
are somehow able to photoionize a significant fraction of 
the ISM not only in the disk but also within the halo, 
1- 2 kpc above the midplane. The heating of the WIM or 
DIG remains poorly understood, however. As shown by 
Madsen, Reynolds & Haffner (2006), the temperatures of 
the Galactic WIM is 2 000 -3 000 K higher than that of 
Galactic Hii regions. Interestingly, the ratios of forbid- 
den lines such as [S ii] , [N ii] or [O i] with respect to Ha 
have been shown to anticorrelate with the WIM/DIG Ha 
surface brightness (e.g. Rand 1998, hereafter R98; Mad- 
sen, Reynolds & Haffner 2006). The fainter the emission 
brightness, the higher these ratios become. A similar 
trend is found with scale height z above the disk in edge- 
on galaxies. Reynolds etal. (1999: RHT) propose that 
an increase in gas temperature may suffice to explain 
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such increase in line ratios. The detailed study by Collins 
& Rand (2001) of the DIG in four edge-on galaxies sug- 
gests, however, that a systematic increase in ionization, 
at least of the fraction of O"*"-!-, might be required as well 
to account for the observed trends. 

Turbulent mixing layers (hereafter TML) induced by a 
hot wind is one of the mechanisms considered by R98 to 
provide supplemental ionization and heating. It allowed 
the author to satisfactorily reproduce the line ratios ob- 
served in the DIG of the edge-on spiral NGC 891 (situ- 
ated at 9.5 Mpc). In order to obtain a reasonable fit, how- 
ever, the mixing layer model had to be combined with a 
varying proportion of matter-bounded photoionized con- 
densations. In this composite model, the photoionization 
models were taken from J. Sokolowski (more information 
can be found in Bland-Hawthorn etal. 2001) while the 
mixing layer calculations were borrowed** from Slavin, 
ShuU & Begelman (1993: hereafter SSB). In the current 
work, we propose two improvements to the heuristic com- 
posite model introduced by R98: a) an integrated model: 
the mixing layer and the matter-bounded photoionized 
condensations are calculated as a single gas component 
submitted to an external ionizing ffux, b) better turbu- 
lence microphysics: rather than considering a single TML 
temperature given by the geometrical mean of the warm 
and hot phases, we derive a temperature structure for 
the TML using the mixing length scheme presented by 
Canto & Raga (1991), later implemented in the emission 
line code mappings ic (Binette etal. 1999). 

The paper is organized as follows. The equations de- 
scribing a turbulent mixing layer and their implementa- 
tion in the code mappings ic are discussed in §[2] The 
results obtained from calculations of the TML structure 
and of its emission line spectrum are presented in §[3] and 
compared with the line ratios from the DIG of NGC 891. 



In particular, the model of SSB with vt- 
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:25 km s 



2 



Binette et al. 



A brief discussion follows in §|4l 

2. MODELING TURBULENT MIXING LAYERS 

2.1. Previous works on TML 

It is not so clear what exactly the difference is between 
"standard" shock wave heating (i.e., heating and/or ion- 
ization due to the passage of a single, well defined shock 
wave), and the heating that takes place in a compressible, 
turbulent mixing layer. In particular, it is not entirely 
clear whether such a turbulent layer actually generates a 
number of weak shock waves, or whether it has a turbu- 
lent cascade ending in turbulent dissipation without the 
appearance of shock waves. 

The initial work on the theory of astrophysical tur- 
bulent mixing layers was presented by Kahn (1980) in 
the context of Herbig-Haro jets. Kahn studied the lin- 
ear and quadratic perturbation theory of the interface 
at the edge of a jet flow. The linear theory was stud- 
ied by a number of authors in considerable detail in the 
context of extragalactic jets in both the gas dynamic 
and the magnetohydrodynamic contexts (see, e.g., the 
review of Bodo 1998). For the radiative Herbig-Haro jet 
case, numerical simulations have been carried out, e.g., 
by Rossi et al. (1997) and by Stone, Xu & Hardee (1997). 
Analytic models based on the standard "turbulent vis- 
cosity" approach (with a turbulent viscosity parameter- 
ized with a simple, "mixing length" approximation) were 
computed by Canto & Raga (1991) and Noriega- Crespo 
et al. (1996). Also, Dyson et al. (1995), Lizano & Giova- 
nardi (1995) and Taylor & Raga (1995) computed mod- 
els that included more or less detailed treatments of the 
chemistry associated with mixing layers. However, these 
models did not include a proper calculation of the ion- 
ization state and line emissivity of the ionized regions 
of these flows. Binette etal. (1999) adopted the mix- 
ing layer scheme developed by Raga & Canto (1997) and 
incorporated the ionization evolution of all the species 
considered and the integration of the full emission line 
spectrum from the mixing layer. These initial calcula- 
tions were intended to explain the spectrum from small 
excitation knots observed in Herbig-Haro objects. 

In the context of of the interstellar and the intraclus- 
ter medium, Begelman and Fabian (1990) proposed a 
prescription for evaluating the temperature of a turbu- 
lent mixing layer. In a subsequent work, SSB used the 
same prescription to calculate the emission line spectrum 
from such layers. Their model included radiative trans- 
fer, which allowed these authors to consider the effect of 
the ionizing photons emitted by the mixing layer itself. 
These calculations were used by R98 in his study of the 
NGC891 DIG. 

For the current work, we adopted the mixing length 
prescription of Binette et al. (1999), described below. We 
also incorporated radiative transfer, which will allow us 
to considers the impact of ionizing photons that originate 
not only from the internal diffuse field produced by the 
mixing layer, but from external sources as well, such as 
the ionizing UV from H II regions. 

2.2. A mixing length approach 

For simplicity, we assume an infinite plane interface, 
along which a hot wind of temperature is flowing 
supersonically with velocity Vw with respect to a static 




Fig. 1. — Schematic diagram showing the cross-section along 
the y-axis of a plane-parallel mixing layer of thickness lo ■ The x- 
axis (parallel to the flow) is structure-less in this representation. 
The mixing layer is formed by the interaction of a supersonic hot 
wind (of temperature and velocity Vw) with a static layer of 
warm gas (of H number density rio and temperature To). The gas 
is entrained into the mixing layer and possess a velocity v, which 
increases monotonically from v = (at j/ = 0) up to Vw The 
different layers are all isobaric, with pressure Po (~ 2knoTo). The 
static layer may contain a photoionized region labeled "Hll em." 
in the figure. Two sources of ionizing photons can be considered: 
external UV sources of flux due to hot stars and the X-ray 
diffuse field generated by the hot wind Lpx (of column A^x and 
average temperature Tx)- The relative importance of the turbulent 
and static layers is a free parameter. 

warm gas layer of hydrogen number density rio and tem- 
perature To- The turbulent mixing layer that devel- 
ops between the hot and warm phases has a geometri- 
cal thickness Zq, as described in the diagram of Fig.[TJ 
For the case of a thin, steady state, high Mach number 
radiative mixing layer, the advective terms along the di- 
rection of the mean flow can be neglected with respect 
to the corresponding terms across the thickness of the 
mixing layer (see Binette et al. 2009 for a more detailed 
description). Under this approximation, the momentum 
and energy equations can be written as: 

where y is a coordinate measured from the onset of the 
mixing layer (see Fig. [1]) , v the bulk flow velocity along 
the x-axis as a function of y, L and G are the radia- 
tive energy loss and gain per unit volume, respectively, 
and /i and k are the turbulent viscosity and conductivity, 
respectively, which are assumed to be constant through- 
out the cross-section of the mixing layer. Eq.[T] can be 
integrated to obtain the linear Couette flow solution: 

v{y) = f K, , (3) 

where lo is the thickness of the TML and Vyj the mean 
velocity of the hot wind (directed parallel to the static 
warm phase's surface). This solution can be substituted 
in Eq.[2l which can then be integrated to derive the tem- 
perature cross-section across the layer. 

Within the TML, we consider the equations govern- 
ing the fractional abundance fi of each species i. This 
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abundance must satisfy the equation: 



D- 



(4) 



where 5*^ is the net sink term (including coUisional ion- 
ization, radiative and dielectronic recombination, charge 
transfer... and processes which populate the current 
species) of the species i. The turbulent diffusivity, D, is 
of order unity and assumed to be position-independent. 
At the inner and outer boundaries of the mixing layer, 
the ionization fractions are set by the equilibrium values. 

To complete the description of the mixing layer, we 
require lateral pressure equilibrium (which determines 
the density of the flow along the axis y), and calcu- 
late the turbulent viscosity with a simple, mixing length 
parametrization of the form: 



(5) 



where pi and ci are the mass density (gcm""^) and sound 
speed (respectively) averaged^ over the cross-section of 
the mixing layer, lo is the geometrical thickness of the 
mixing layer (Fig. [1]). 

From the work of Canto & Raga (1991), the value of 
the proportionality constant a is 0.00247, as clarified in 
Appendix]^ This is the required value for a supersonic 
mixing layer model to match the opening angle of « 11° 
of subsonic, high Reynolds number laboratory mixing 
layers in the limit in which the jet Mach number tends 
to one. Canto & Raga (1991) also find that mixing layer 
models with this value for a reproduce the opening angle 
of mixing layers in jets with Mach numbers^, M^,, of up 
to 20. Our treatment of the TML can be extended to 
cover the subsonic case by replacing ci in Eq.[5] by Vw , 
whenever Vw < ci (see AppendixlA]). 

Considering that the turbulent conduction and diffu- 
sion Prandtl numbers are of order one, we can compute 
the conduction coefhcient (where Cp is the 

heat capacity per unit mass averaged across the mix- 
ing layer cross-section) and the diffusion coefficient as 
D Ki n/p. In this way, we obtain a closed set of second 
order differential equations, Eqs. ([2]) and (jl]), which can 
be integrated with a simple, successive overrelaxation nu- 
merical scheme. 

2.3. TML calculations with the multipurpose code 
MAPPINGS IC 

We use the code mappings ic (Ferruit etal. 1997) 
to compute the radiative energy loss term L and the 
photoheating term G (Eq.[2|) at each position across the 
TML. At both the inner and outer boundaries of the 
mixing layer, we assume equilibrium ionization of the 
different species, while across the layer, our simple over- 
relaxation scheme is used to determine the ionization 
fractions (Eq.|3]). For the ion diffusion of each species 
fi, the spatial differential equations are converted to 
temporal equations, with the use of pseudo-time steps 
At = Ay^ lo/{cy.nci), where n is the average H number 
density. This allows us to use the temporal algorithm 

^ The subindex I denotes averages across the mixing layer, see 
AppendixIXl 

It is customary to express V-w in terms of the wind Mach num- 
ber, Mw = Vw/co, where Co is the adiabatic sound speed in the 
static layer. 



previously described in Binette & Robinson (1987) for 
determining the spatial diffusion of the ionic species. 

The radiative transfer is determined by integrating 
(from the hot layer y = lo down to y = Q) the intensity 
of the UV diffuse field that the layer produces, assuming 
the outward only approximation. Any external UV radi- 
ation impinging from the external side is simply added 
to the diffuse field at the onset of the integration. The 
intensity of the external ionizing field is defined by the 
ionization parameter as follows: 



(6) 



where c is the speed of light, the fiux number of ioniz- 
ing photons impinging on the TML, and rio, the total H 
density at y = 0. The flux ip^w is the sum of two compo- 
nents: resulting from external UV sources such as hot 
stars and ip>x due to the X-ray diffuse field generated by 
the hot wind (of column A^x and average temperature 
Tx). The absorption processes that are taken into ac- 
count in the transfer equations across the TML include 
all the relevant photoelectric cross-sections of the ions 
present. 

Collins & Rand (2001) used a method based on 
[Nii]A6583 to estimate the gas temperature and, by boot- 
strapping, they could infer the ionization fractions of N 
and other species. They find that the hydrogen ioniza- 
tion fraction in NGC891 is as high^ as 0.55-0.80. Be- 
cause H is fairly ionized and since X-ray ionization is 
very inefficient in comparison with UV radiation, we in- 
fer that their role relative to the soft UV fiux from stars 
must be relatively unimportant. We therefore neglected 
the term ipx by setting A'x to a negligible value. 

For all the calculations presented in this work, the mix- 
ing layers are considered immersed in the radiation field 
that escapes from Hii regions located in the disk. This 
picture assumes that disk H ii regions are on the whole 
matter-bounded. According to Beckman etal. (2000), 
Zuritaetal. (2002) and Relafio et al. (2002), the fraction 
of Lyman photons that escape giant H ii regions (labeled 
"leakage fraction" in ij |3.2p may lie in the range of 30- 
50%. The impact of photoionization on our TML model 
can be inferred from the internal behavior of F across the 
layer, where F is defined as follows: 



F = 



L-G 
L + G 



(7) 



The quantity F is zero when the temperature corresponds 
to the equilibrium value and tends toward unity when 
cooling becomes much stronger than heating. In our 
calculations, we iteratively determine To at y = un- 
til F = is obtained. 

After computing the emission line spectrum of a given 
TML, MAPPINGS IC offers the option of separately com- 
puting the emission lines generated by the inner pho- 
toionized layer (i.e. y < in Fig.[l]) where equilibrium 
ionization prevails. A simple isobaric photoionization 
model is calculated in this case, using the radiation field 
that has not been absorbed by the mixing layer. The 
total line spectrum from all the layers is then given by 
the sum of the line intensities from both the TML and 
the static equilibrium photoionization calculation. For 

Comparable values were inferred for the DIG in the Milky Way 
by Reynolds etal. (1998) using [Ol]A6300. 
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all the models studied, however, whenever a static pho- 
toionized layer was included, it dominated the total spec- 
trum and dwarfed the spectral signature of the turbulent 
layers. In order to enhance the contribution of the TML 
and successfully fit the ratios observed in NGC891, the 
static equilibrium region must be obviated and the to- 
tal thickness of the emission layer becomes equal to that 
given by the turbulent layer alone. This might come 
about in a 3D geometry if the cold gas were to reside 
in the core of dense condensations whose area is much 
smaller than that of the surrounding turbulent layers. 

Since the mixing layer is isobaric, with pressure Po, 
the density profile depends on the behavior of the tem- 
perature within the layer itself. The value of Po is set 
by two quantities evaluated at the inner TML boundary 
{y — 0): the density Uo and the equilibrium tempera- 
ture To- To be definite, a density of rio = 1 was adopted 
that applies to all models presented in this work. At 
such a low density, the low density regime fully applies 
and the calculated line ratios are therefore independent 
of the actual rio value used, provided the product nolo is 
conserved. Our models are appropriate for any density 
^ 50. When the external parameters {Uo, T„ and V^) 
are kept constant, TML models turn out to be equivalent 
whenever the product of the H density rio and thickness 
lo is the same. Therefore, at least in the low density limit 
regime typical of the DIG, it is sufficient to specify the 
quantity® No=nolo to uniquely define a model. 

We found that the temperature of the hot wind is not 
a critical parameter of the TML models. Its value can be 
varied and still result in a sequence of equivalent models 
(i.e. with similar line ratios), provided the value of 
is kept constant. For this reason, we adopted the same 
value of = 8 X lO'^ K for all the models presented in 
this work. 

To summarize, in order to compute solutions to the 
mixing layer, we must first define the external radiation 
field and then specify the values of the following param- 
eters: the mixing layer's No, the ionization parameter, 
Uo, and finally the velocity of the hot wind, V^, and its 
temperature, which we kept constant at 8 x 10^ K. 

3. MODELING THE DIG IN NGC 891 

3.1. The gas abundances 

Unless otherwise specified, the abundances of the el- 
ements are solar (Anders & Grevesse 1989) and have 
the following values: (He, C, N, O, Ne, Mg, Si, S, Ca, 
Ar, Fe) = [10^, 363, 112, 851, 123, 38, 36, 16, 2.3, 3.6, 47] x 
10~^ relative to H (by number). Depleted or oversolar 
metallicities resulted in line ratios difficult to reconcile 
with the forbidden line ratios observed in NGC 891. 

3.2. The ionizing spectral energy distribution 

The first step consists of defining the spectral energy 
distribution (sed) of the ionizing radiation field reach- 
ing the DIG. A useful constraint is provided by the 
observed He i A5876/IIa ratio. In a previous observa- 
tional work on NGC 891, Rand (1997) measured a ratio 
of HeiA5876/Ha « 0.034, which he estimated implies 

* The quantity A^o is a convenient model descriptor. However, 
it is a bad estimator of the true integrated H column, which is a 
lot smaller since the density is not constant but decreases as the 
temperature rises with thickness y (Fig.[T]l. 
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Fig. 2. — The sed of a 38 000 K star (dashed hue) as a function 
of photon energy. The solid line represents the SED of the UV 
radiation escaping from a matter-bounded photoionized slab that 
absorbs 80% of the impinging ionizing photons. It is assumed that 
the DIG is exposed to this absorbed SED. 

a ratio of He ionizing photons to H ionizing photons, 
Qh/Qho, of about 0.08. To be consistent with this da- 
tum. Rand favors a stellar temperature of 37 500 K. In 
the calculations presented here, we adopt a similar tem- 
perature of 38 000 K. The atmosphere models used in 
MAPPINGS ic are from Hummer & Mihalas (1970). As in 
R98, we assume that only a fraction of the stellar radia- 
tion escapes the H ii regions in the spiral disk to reach the 
DIG. To be definite, we have set the leakage fraction to 
20%. To compute this modified SED, we have computed 
a plane-parallel photoionization model and extracted the 
photon spectrum at a depth where 80% of the ionizing 
photons have been used up. The diffuse ionizing radia- 
tion from the slab itself has been included. The resulting 
SED is shown in Fig. [2] (solid line). Because this modified 
SED contains few ionizing photons beyond 24.6 eV, the 
mixing layer model presented below is characterized by 
an He i A5876/Ha ratio as low as ~ 0.003. This suggests a 
higher stellar temperature for NGC 891 than considered 
here and in Rand (1997), or a larger leakage fraction than 
assumed above. 

We have verified that the stellar temperature can be 
raised up to 40 000 K without exceeding the observed 
Hci/Ha value (assuming the same leakage fraction). For 
the DIG in the Milky Way, HeiA5876/Ha is observed 
to be significantly smaller, lying in the range 0.01-0.04 
(Madsen etal. 2006). 

3.3. TML temperature and ionization structure 

In Fig. Owe present the behavior of the temperature as 
well as of the imbalance between cooling and heating (F, 
see Eq.[7]) as a function of thickness y within a represen- 
tative mixing layer calculation. For this model, the hot 
wind velocity is 84kms~^, while the ionization parame- 
ter is Uo = 10^^ and the layer's column iVo — lO^^cm^^. 
The inner boundary is characterized by an equilibrium 
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Fig. 3. — Behavior of T (Eq.[7ll and of the local temperature as a 
function of (a) the electron column density (normalized to the total 
value of 1.1 X 10^® cm~^) and of (b) the normalized thickness y/lo- 
The model shown has y„=85kms-i, T„ = 8 X 10^ K, Uo = 10"^ 
and No = lO^'^cm"-^. The temperature scale is normalized relative 
to Tw The value of the equilibrium temperature, at the onset of 
the layer, is ro=4390K (i.e. To/Tm = 0.0055). 
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Fig. 4. — Ionization fractions of O, S, He and H for the same 
model as in Fig. [3] as a function of (a) the electron column density 
(normalized to the total value of 1.1 X 10^® cm"'^) and of (6) the 
normalized thickness y/lo- Also shown is the quantity (n X rie)^ 
(thick grey dashed line), which is normalized to the value at j; = 
and is plotted on a logarithmic scale, to be read using the rightmost 
axis. The volume line emissivities from each ion is proportional to 
this factor. The lines' coding is the same in both panels. 



temperature of To = 4 390K. As illustrated in Fig.O due 
to the very hot wind temperature (T^, » To), T increases 
monotonically to reach unity at y — lo &s a result of the 
cooling term becoming progressively dominant relative 
to the photoheating term. 

The ionization structure within the same TML is 
shown in Fig. HI The static layer if it was present would 
have ionization fractions quite similar to those found at 
the extreme left in this figure. Line intensities depend on 
temperature, but scales also with density square. Since 
the quantity (n x ne)^ decreases strongly with increasing 
thickness (see thick grey dashed line), the local emissiv- 
ities likewise decrease markedly towards the right. It is 
therefore apparent that any emission from [O ill] or any 
other species with similar or higher ionization potential 
are produced in this model only within the mixing layer, 
and not within the static layer. This would of course be 
different if a significantly higher stellar temperature was 
chosen. 

3.4. DIG line ratio trends in NGC891 

For comparison with TML models we selected the work 
of R98 on the DIG in NGC 891. The data on this object 
present the following advantages: the spectral data are 
homogeneous, they are of sufficient quality to reveal sub- 
tle intrinsic correlations among the line ratios observed 
at different positions across and along the disk of this 
edge-on spiral, finally, mid-infrared halo and disk obser- 
vations with the Spitzer Space Telescope are now avail- 
able (Rand etal. 2008). We postulate that the line ratio 
correlations are a direct manifestation of the DIG exci- 
tation mechanism rather than positional changes in gas 
metallicity or in effective stellar temperature. 

To optimize the comparison with models, we present 
the line ratios in Fig. [5] in a format which reveals 
the essential trends found by R98 in his long slit 
data. The three ratios [Oiii]/H/3 (5007/4861), [Oi]/Ha 
(6300/6563) and [Nii]/[Sii] (6583/6716) are aU plot- 
ted against the same ratio [Sll]/Ha (6716/6563). The 
[Oiii]/H/3 line ratio lies below unity in all measurements 
(panel a) . What is most striking about the data is that 
not only do the forbidden line intensities [Oi], [Nil] 
increase along with [Sii]/Ha, but the high excitation 
[Oili]/H/3 ratio as well. Such a behavior runs counter to 
the predictions of simple photoionization models where 
[Oiii]/H/3 increases markedly when the ionization pa- 
rameter is increased, while the ratios [Sii]/IIa, [Oi]/Hq; 
or [N ii]/IIa all decrease. This leads to the proposition by 
R98 and Collins & Rand (2001) that another ionization 
process besides photoionization is likely to be involved. 
With the Milky Way DIGs, the derived temperatures are 
found to increase along with the increase of the forbid- 
den lines relative to Ha or II/3 (e.g. Madsen etal. 2006). 
According to RHT, a systematic temperature increase 
might be all that is necessary to reproduce the increase in 
the forbidden line intensities. Again, this runs counter to 
predictions of simple photoionization models where the 
temperature is quite insensitive to changes in the value 
of Uo- Hence various authors have concluded in the exis- 
tence of a new heating (e.g. RHT) and possibly ionizing 
mechanism (e.g. R98) in addition to photoionization. 
For the sake of completeness, we included an inset in 
Fig. [5] (panel d) where the [N ii] temperature computed 
in models is plotted against [Sii]/Ha. We may reason- 
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Fig. 5.— Line ratios of (a) [Oni]/H/3, (b) [Oi]/Ha, (c) [Nii]/[Sii] 
and (d) the Tjmh temperature, all as a function of the [Sll]/Ha 
line ratio. Open squares: measurements of R98 of the DIG in 
NGC891; continuous line: a sequence of TML models in which 
Vw is progressively increased by 0.1 dex from 27 to 133kms~^, 
while other parameters are kept constant: No = 10^^ cm~^ and 
Uo = 10~^. The ratio [Sll]/H« increases with increasing Vw The 
dotted line represents a varying mixture of turbulent mixing layers 
and photoionized matter-bounded models that were proposed by 
R98. 

ably expect that any realistic DIG models should show 
a significant temperature gradient in such a diagram. 

3.5. Specific TML calculations 

We carried out an extensive exploration of the param- 
eter space in order to determine the parameter range 
that could best reproduce the trends present in Fig. [51 
As in previous studies, we concur that the layer has to 
be matter-bounded, otherwise the [O i]A6300/[S ii]A6716 
ratio would significantly exceed the observed values. Fur- 
thermore, for the calculations to show a significant de- 
parture from steady-state photoionization as a result of 
turbulent dissipation and ionization mixing, it required 
an extremely thin layer, with a column iVo of order 
~ 10^^ cm~^ only. This raises concerns about the effi- 
ciency of turbulent dissipation to contribute significantly 
to the DIG emission, as discussed in §|31 Because the 
TML layer in the composite model of R98 was decou- 
pled from the photoionized component, the author was 
free to adopt a much larger column of 2 x 10^^ cm^^ in 
his matter-bounded photoionization calculations. 

A characteristic of the DIG is its relatively weak [O ill] 
line intensity, which favors relatively small values for the 
ionization parameter when considering photoionization 
(e.g. Sokolowski & Bland-Hawthorn 1991). In the case 
of the composite TML models of R98, the author has 
employed values that varied between 10^^ to 8 x 10^^. 
In the TML regime that we have explored, we are lead 
towards similar values, with Uo typically at ~ 10~^. We 
find that the [Oiii] flux is the result of turbulent dis- 
sipation rather than of photoionization. On the other 
hand, low ionization species responsible for lines such as 



Fig. 6. — Dashed line: a sequence of TML models in which the 
ionization parameter is progressively increased by 0.2 dex, from 
4.0 X 10-*^ to 4.0 X 10-5. The column A^o increases in the same 
proportion, from 0.4 X lO^'^ to 4.0 X lO^'^ cm" ^ while Vw is kept con- 
stant at SSkms" The ratio [Sll]/Ha decreases with increasing 
Uo and A^o. Other symbols have the same meaning as in Fig. [5] 

[On], [Nil], [Sii]... are mostly caused by photoioniza- 
tion, but the additional heating due to turbulent dissi- 
pation causes these lines to become much brighter. In 
fact, if a static layer (see Fig.[T]) with such a low value 
of Uo was present behind the TML, the overall line ra- 
tios would lie near [N ii]/[S ii]=0.93, [S ii]/Ha=0.087 and 
[O i]/Hq:=0.013. Therefore, these ratios would all lie out- 
side the boundaries from the three line diagrams a, b and 
c of Fig.O If we relied only only photoionization (with- 
out any TML) to fit the observations, we would have 
to adopt a much harder ionizing SED (but see ij |3.2p to 
provide the necessary heating and the ionization up to 

0+2. 

3.5.1. The Vu,- sequence 

In Fig.[5l we superpose a sequence of TML models to 
the data, in which the wind velocity is progressively in- 
creased in locked steps, by 0.1 dex, from Vu, = 27 to 
133kms~^ (corresponding to Mach numbers increasing 
from Mtu=2.53 to 13.5). The column Nq was set to 
lO^^cm-2 and Uo to 10~^. In that sequence, the forbid- 
den line ratios of [0iii]/H/3 and [Sii]/H/3 increase with 
increasing Vw We note that the temperature inferred 
from the [Nil] (5755/6583) line ratio, Tnii, increases 
slightly along with [S ii]/Ha, as shown in the upper right 
inset (panel d) . Overall, the trends in line ratios are well 
reproduced, better in fact than with the composite (pho- 
toionization-f TML) model of R98, which is represented 
by the dotted line. R98 showed that the line ratios not 
only correlate with decreasing Ha surface brightness, but 
also with increasing z height above the disk of NGC 891. 
If the Via sequence of models were a valid description of 
the long-slit measurements of R98, this would imply that 
the wind velocity is increasing with z. 
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-0.4 -0.2 
log [NIl]/Ha (6583/6563) 

Fig. 7. — The Tnii temperature was derived from the temper- 
ature sensitive [Nil] (5755/6583) line ratio as a function of the 
[Nil] /Ha line ratio. The three solid squares correspond to temper- 
atures borrowed from of Madsen et al. (2006) and Reynolds et al. 
(2001) of the DIG in the Milky Way while the two upper limits 
correspond to NGC891 (Rand 1997). Overlayed are the same Uo 
and Vw sequences of Fig. [5] (solid line) and Fig.|6] (dashed-line) , 
respectively. 

3.5.2. The Uq- sequence 

Alternatively, it is possible to reproduce the line ratio 
correlations by varying the ionization parameter instead 
of • In Fig-IH we show a sequence of models in which 
both Uo and A'o increase^ in lock steps by 0.2 dex, while 
Vw is kept constant at 85kms~^. The initial values of 
Uo and iVo are 4.0 x 10~^ and 0.4 x lO^'^cm"^, respec- 
tively, while the final values are 10 times larger. Varying 
No in locked steps with the ionization parameter limits 
the variations in H i opacity, which is desirable otherwise 
the TJVIL layers either get too thin or reach a regime in 
which photoionization takes over. In the sequence shown 
in Fig.[6l we have ensured that photoionization does not 
become the dominant process. In fact, [Oiii] and the 
other forbidden lines actually decrease as Uo is increased. 
For instance, [O ill] gets weaker as Uo increases because 
the average gas temperature where it is produced de- 
creases along the sequence from 51500 to 40 300 K. This 
sequence of models shows a [Nil] temperature gradient 
(Fig.[6H) that is stronger than in the previous sequence 
(Fig.EH). When compared to the data on NGC891, the 
behavior of the [O ill] /II/3 ratio is not as closely fit as in 
the Ku- sequence, but the fit of the other line ratios are 
better. 

3.6. Other DIG observations 

Although this study focuses on NGC891, fitting simi- 
lar line ratio trends in other objects is certainly achiev- 
able. In itself, this does not validate TJVIL models, since 
the main reason for this is the higher number of free 

^ Such a sequence is achieved by increasing both quantities ipw 
and lo in lock steps, while keeping Uo constant. 



parameters that can be varied and which were already 
numerous in the alternative case of pure photoionization. 

Although forbidden line ratios with respect to Ha 
or II/3 increase with height z above the midplane in 
NGC891, this behavior is not repeated in all spirals. 
Cases for instance where [O ill] /II/3 decreases or reverse 
trend with height z are well documented in the studies of 
Miller & Veilleux (2003), Otte etal. (2002) and Collins 
& Rand (2001). In their studies of face-on M33, Voges 
& Walterbos (2006) find sudden changes in this ratio. 
Within the perspective of TML models, these observa- 
tions bring new constraints to the interpretations that 
models are able to provide. For instance, in cases where 
[O III] /HZ? decreases with z, we find it improbable that Vw 
could be decreasing with height, as required by the V^- 
sequence (Fig. [5]). In the case of the Uo~ sequence (Fig. El), 
however, a TML interpretation would imply that Uo is 
increasing with height, which could only come about if 
the ambient gas pressure decreases faster with z than the 
local ionizing flux. 

The direct measurement of the [N ii] temperature is dif- 
ficult owing to the intrinsic weakness of the [Nii]A5755 
lines. The [Nil] (5755/6583) line ratio has nevertheless 
been determined in the Milky Way for three DIG sight 
lines (Madsen etal. 2006; Reynolds etal. 2001). The 
inferred temperatures are plotted in Fig.[7l The two up- 
per limits for NGC891, as determined by Rand (1997), 
are also plotted. The comparison of these three temper- 
atures by Madsen et al. (2006) with those of 'classical' 
H II regions confirms that the DIGs are approximately 
2000- 3000 K warmer than Hii regions. The TML mod- 
els presented in this work are overlayed in Fig.[7l Despite 
the scarce data, the comparison with models is encour- 
aging insofar as it shows that turbulent dissipation and 
mixing can provide the extra heating required by the 
Galactic DIGs. 

With the TML models presented in this work, in which 
none of the [O ill] emission is the result of photoioniza- 
tion but rather due to turbulent mixing and dissipation, 
we find that [O ill] is produced at a depth y of very 
high electron temperature. The behavior of the [O ill] 
(4363/5007) line ratio in models shows that in both se- 
quences the increase in [O ill] is caused by an increase in 
average^" temperature, from ~ 40 000 to 52 000 K. These 
values are much higher than those characterizing the 
[Nil] emission zone, to the extent that the [Oiii]A4363 
line flux is between the value of the [Nii]A5755 flux and 
half of it. In principle, measuring a high [Oiii] temper- 
ature or not would provide a strong discriminant for or 
against TML models. 

3.7. The mid-infrared [Nell] 12.81 fj,m and 
[Neiii] 15.56 fim lines 

Mid- infrared observations of NGC891 using the 
Spitzer Space Telescope has been reported by Rand et al. 
(2008). Of particular interest are the measurements of 
the [Neil] 12.81 ^m, [Neiii] 15.56 /.tm and [Siii] 18.71 /.tm 
fluxes of the DIG in the halo (2 measurements) as well as 
in the disk (1 measurement). Interestingly, the observed 

Because of the steep temperature gradient across the TML, 
the integrated [Olll] (4363/5007) line ratio represents an 'average' 
temperature that is strongly biased towards layers where the gas 
is densest rather than being close to the temperature where the 
fraction of O^"^ peaks (see behavior of (n X rie )2 in Fig.lUl. 
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[Ne III] /[Neil] ratio increases with height z by a factor of 
three. This is quahtatively similar to the trends revealed 
by the optical ratios of [Sii]/Hq; or [Oiii]/H/3. Our two 
TML sequences has [Neiii]/[Neii] increasing along with 
[Oiii]/H/3, but the increase is only a factor two. The 
maximum values reached for the ratio [Neiii]/[Neii] are 
0.067 and 0.070, for the 14,- sequence and C/q- sequence, 
respectively. This is about 4.5 times too weak with re- 
spect to the observed halo value of 0.31. In his analysis. 
Rand etal. (2008) adopted a Ne/0 abundance 2.5 times 
higher ratio than the solar value (ij lB.ip . Changing our 
abundance set might resolve this discrepancy. 

4. DISCUSSION 

A photoionized turbulent mixing layer can reproduce 
the line ratio trends observed by R98, if it is matter- 
bounded. Our work confirms the initial intuition of R98 
who mixed two causally unrelated ionized components 
in ad hoc proportions: a TML component and a matter- 
bounded photoionized component. The main advantage 
of our integrated model is that the ionized slab consid- 
ered is submitted in a self-consistent manner to both 
excitation mechanisms: turbulent dissipation and pho- 
toionization. Adjusting the parameters so that the ob- 
served line ratios can be reproduced, however, comes at 
a price. The thickness of the TML layer becomes so 
small that the overall efficiency in generating emission 
lines becomes a concern. This can be illustrated as fol- 
lows. From the definition of Uo, the ionizing photon 
flux impinging on our putative mixing layer of area A 
is 3.0 X 10'^ An^i U-z quantas^\ where [/_5 is f/o/lO^^ 
and n^i = rio/O.l cm~^. In the case of pure photoioniza- 
tion, assuming ionization equilibrium, the Ha luminosity 
is given by Lhq = 4.44 x lQ~^'q Ari-i U-^evgs~^ where 
?7 is a measure of the reprocessing efficiency, that reaches 
unity in the case of an ionization-bounded layer, where all 
the ionizing photons are absorbed. When the Lhq from 
the TML models are compared to that expected from 
this simple relationship, we find that the 14;- sequence 
is characterized by an efficiency rj that goes from 0.013 
to 0.0031, as Vw is increased. In the case of the Uq- 
sequence, 77 goes from 0.0013 to 0.027, as Uo and No 
are increased. What is most problematic, is that the effi- 
ciency is lowest when turbulent dissipation is maximized, 
that is, in the regime where the [Oiii]/H/3 ratio is high- 
est, whether we consider the 14,- sequence of Figs. [5] or 
the C/q- sequence of Fig. HI The efficiencies at the high 
[Oiii]/H/3 end are only 0.0031 and 0.0013, for the V^- 
sequence and ?7o- sequence, respectively. 

Such an efficiency is uncomfortably low. This can 

For this particular model, the layer's geometrical thickness is 
cm, which in any case leaves no space for the 



be illustrated by using the observations of NGC891 
by Rand (1997). With the long-slit perpendicular 
to the galactic plane of edge-on NGC891, the au- 
thor measured an Ha surface brightness of 2 x 10^^^ 
erg cm~^ s~^ arcsec"^ at an height 2;=2kpc, where the 
[Oiii]/H/3/ ratio peaks. Assuming a distance oi D = 

9.5 Mpc (1 arcsec=46pc), the slit, perpendicular to the 
galactic plane, crosses the mid-plane of NGC891 at 

4.6 kpc from the nucleus. We adopt 1" as the refer- 
ence spatial scale. We integrate along the thickness of 
the layer, which we assume extends at most over 9 kpc 
along the line of sight, so that the layer area A' ex- 
posed to the UV flux arising from the midplane be- 
comes 46pcarcsec^^ x 9 000pc. The surface brightness 
derived from the Lfja expression above becomes: Fhq = 

k A' / (AttD^), where k is the number of projected 
TML layers per arcsec (at z—2 kpc) along the slit. Hence, 
Fhq = 1.6 X 10~^^fc 7yn_i t7_5 erg cm~^ s~^ arcsec"^. For 
an efficiency rj of 0.0013 and with C7_5=0.4, matching 
the observed surface brightness implies the unreasonable 
number of fc ~ 2400 layers^^ per arcsec along the slit, 
assuming a density typical of the DIG in the Milky Way 
(n_i=l). Even if one raises this density by an order of 
magnitude, too many layers are still required. For higher 
density values beyond that, there would not be enough 
O stars to account for the required ionizing flux. 

To summarize, it is possible to contrive the values of 
the parameters describing TML models so that the ob- 
served trends in line ratios are reproduced reasonably 
well. However, the efficiencies of any successful models 
appear downright too insignificant. In its current form 
our TML model is therefore invalidated. We may spec- 
ulate that the gradient in temperature between the hot 
wind and the warm photoionized gas is so steep that the 
mixing length approach explored by Canto & Raga is not 
directly applicable to the very hot wind surrounding the 
DIG. Values for the constant a (Eq.[5|) about two orders 
of magnitude larger would be required so that turbulent 
mixing and dissipation could operate over a layer's thick- 
ness much larger than the very small values found in the 
current model. 
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hot wind to flow if k is large. 



APPENDIX 
THE DERIVATION OF CONSTANT a 
From Eq. A3 in Canto & Raga (1991), we obtain that the turbulent viscosity coefficient fi is given by: 

fJ-^ Y^PoColo , (Al) 

where e is given by 0.089 x MIN[ei, €2] with ei = and £2 = (see their Eqs. 27, 30 and 36). c^, and Co correspond 
to the adiabatic sound speed in the hot wind and in the static layer, respectively. The subindex I denotes averages 
across the mixing layer. Canto & Raga (1991) refer to the £2 < £1 case as the 'mixing layer limited' case while the 'jet 
limited' case corresponds to £2 > £1 
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In the case where £2 < ci, since in the isobaric case picf = PoCg, after multiplying the numerator and denominator 
by pici, we can rewrite the above expression as 

0-089 , 

In the notation used in this Paper, the constant a simply corresponds to the ratio 0.089/36 = 0.00247. In the transonic 
case or when ei < £2, we replace q in Eq. lA2l bv MIN[T4> , Q x MIN{1, |^}]. It turns out that the value of a used by 
Binette etal. (1999) is a factor 3 too large. This can be corrected by dividing by 3 the values of the layer's thickness, 
h, quoted in their paper. 
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